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We investigate the cosmic age problem associated with 9 extremely old globular clusters in M31 
galaxy and 1 very old high-2 quasar APM 08279 + 5255 at z = 3.91. These 9 globular clusters 
have not been used to study the cosmic age problem in the previous literature. By evaluating the 
^^ ' age of the universe in the ACDM model with the observational constraints from the SNIa, the 

Cn I BAO, the CMB, and the independent Ho measurements, we find that the existence of 5 globular 

^ . clusters and 1 high-z quasar are in tension (over 2a confidence level) with the current cosmological 

Q ' observations. So if the age estimates of these objects are correct, the cosmic age puzzle still remains 

^T' I in the standard cosmology. Moreover, we extend our investigations to the cases of the interacting 

^^ , dark energy models. It is found that although the introduction of the interaction between dark 

CO ' sectors can give a larger cosmic age, the interacting dark energy models still have difficulty to pass 

the cosmic age test. 
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I Before the great discovery that our universe is undergoing an accelerated expansion [l|, |2], the most popular 

O ■ cosmological model was the SCDM model (i.e., a flat universe with the present fractional matter density flm — !)■ 
■*— > . However, the SCDM model is always plagued by a longstanding puzzle: in this model the present age of the universe is 
C^ ^0 — -^jj- — 9 Gyr {Hq denotes the Hubble constant), while astronomers have already discovered that many objects are 
older than 10 Gyr Q. This so-called cosmic age problem becomes more acute if one considers the age of the universe 
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I. INTRODUCTION 



at a high redshift. For instance, a 3.5 Gyr-old galaxy 53W091 at redshift z — 1.55 and a 4 Gyr-old galaxy 53W069 at 
z = 1.43 are more difficult to accommodate in the SCDM model [4]. Along with the discovery of accelerated expansion 
\jr\ ' of our universe and the return of the cosmological constant A, the age problem has been greatly alleviated. The 7- year 
■Tilj- I WMAP observations Q tell us that in the standard cosmological model (i.e. the ACDM model) the present cosmic 
C^ ■ age is to = 13.75±0.11 Gyr. Besides, it is shown that the ACDM model can also easily accommodate galaxies 53W091 
■^ ' and 53W069 '.(y^. Therefore, it is widely believed that the cosmic age problem is a "smoking-gun" evidence for the 
l/~j . standard cosmological model, and many people believe that the cosmic age puzzle does not exist in the ACDM model. 
f^ ■ However, we still want to raise the question: Is the cosmic age problem really removed by the ACDM model? 
^D , In a recent paper [7|, by comparing photometric data acquired from the Beijing- Arizona- Taiwan-Connecticut 

(BATC) system with up-to-date theoretical synthesis models, the ages of 35 globular clusters (GCs) in M31 galaxy 
were estimated by Ma et al. Soon after, the ages of other 104 GCs in M31 galaxy were also obtained by the same 
research group Q. From their results, we find that 9 extremely old GCs are older than the present cosmic age 
(to = 13.75 Gyr) predicted by the 7-year WMAP observations (see Table [T]). To our present knowledge, these 9 
jrt ■ GCs have not been used to study the cosmic age problem in the previous literature; so it is interesting to explore 
■ " implications of these 9 old GCs for the standard cosmology. In addition, the existence of high- 2; quasar APM 08279 

-I- 5255 at z = 3.91 @, [itj] still remains a mystery |10l - [l4j . Using the maximum likelihood values of the 7- year WMAP 
observations flm = 0.272 and h — 0.704 {h is the reduced Hubble parameter) |5|, it is seen that the ACDM model can 
only give a cosmic age t = 1.63 Gyr at z = 3.91, while the lower limit of the quasar's age is 1.8 Gyr [l3|. In a previous 



paper 1J| , Wang and Zhang demonstrated that by simply introducing the dark energy (DE) alone we cannot remove 
the high-z age problem, and suggested that the introduction of interaction between dark sectors may be helpful to 
alleviate the cosmic age problem. But in [lj|, when evaluating the age of the universe in the interacting DE models, 
the parameters describing interaction strength were chosen arbitrarily. A consistent quantitative analysis should use 
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TABLE I: The ages of 9 extremely old Globular Clusters in M31 Galaxy. Notice that the present cosmic age predicted by the 
7-year WMAP observations is to = 13.75 ±0.11 Gyr. 



Object 


Age(Gyr) 


Reference 


B024 


15.25 ±0.75 


Table 5 of [8| 


B050 


16.00 ± 0.30 


Table 5 of [7| 


B129 


15.10 ±0.70 


Table 5 of [7] 


B144D 


14.36 ± 0.95 


Table 5 of [8] 


B239 


14.50 ±2.05 


Table 5 of [T] 


B260 


14.30 ±0.50 


Table 5 of [8] 


B297D 


15.18 ±0.85 


Table 5 of [8] 


B383 


13.99 ±1.05 


Table 5 of [8J 


B495 


14.54 ±0.55 


Table 5 of [8| 



the maxiniuin likelihood method to determine the corresponding model parameters from the current cosmological ob- 
servations. Then using the obtained results one can analyze the consistency with the cosmic age data. This analysis 
will be done in the present work. 

This paper is organized as follows. In section II, we evaluate the age of the universe in the ACDM model with 
constraints from the Type la supernovae (SNIa), the Baryon Acoustic Oscillations (BAO), the Cosmic Microwave 
Background (CMB), and the independent Hq measurements, then use the obtained results to analyze the consistency 
with the age data listed above. In section III, we extend our investigation to the cases of three interacting dark energy 
models. We summarize in Section IV. Appendix A supplies the details of how to utilize the cosmological observations 
to constrain DE models. We assume today's scale factor ao = 1, so the redshift z = a^^ — 1, the subscript "0" always 
indicates the present value of the corresponding quantity. 



II. COSMIC AGE PROBLEM IN THE ACDM MODEL 



The age t{z) of a flat universe at redshift z is [llj 

t{z) = 



dz 



ii + ~z)Hisy 

where H{z) is the Hubble parameter, given by the Friedman equation, 

3MpiH^ = p,n + Pde + Pr, 



(1) 



(2) 



Mpi = I/vSttG is the reduced Planck mass, pm, Pde, and pr are the energy density of matter, dark energy, and 
radiation, respectively. The reduced Hubble parameter is E{z) = H{z)/ Hq. Once E{z) is given, the cosmic age t{z) 
at any redshift z can be evaluated. For ACDM model. 



E{z) = v/^^(l + Z)4 ± n^{l ± z)3 ± 1 - ^m - ^r, 

where fi^ is the present fractional radiation density given by p , 

n^ = n^{l + 0.22717Ve//), n^ = 2.469 X 10"^/^~^ Neff = 3.04. 



(3) 



(4) 



Note that in a flat universe the present fractional DE density is il.de = 1 ^ ^m — ^r • Since we also consider effects of 
radiation component, there are two model parameters, rim and h, need to be determined. 

To constrain DE models by using cosmological observations, we will employ the x^ statistics [la|- For a physical 
quantity ^ with an experimentally measured value ^ots, a standard deviation a^, and a theoretically predicted value 
^th, the x^ value is given by 



Xi 



2 _ {^th — ^obs) 



(5) 
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FIG. 1: Cosmic age problem at z = in tlie ACDM model. The dotted elliptical line, the dashed elliptical line, and the central 
star symbol denote the la confidence region, the 2a confidence region, and the best-fit point of the ACDM model in the Q,m-h 
plane, respectively. Three solid diagonal lines, from top to bottom, represent the constraints of 13.82 Gyr, 13.88 Gyr, and 14 
Gyr at z = 0, respectively. Notice that the top diagonal line tangents to the la confidence region of the ACDM model, while 
the middle diagonal line tangents to the 2a confidence region of the ACDM model. For a given diagonal line, the area below 
this diagonal line corresponds to a larger cosmic age. The ACDM model predicts a lower cosmic age than the la lower limits 
of 5 GCs' age. 



The total x^ is the sum of all x?s, i.e. 






(6) 



The model parameters yielding a minimal x^ is favored by observations. The observational data used in this paper 
include the Constitution SNIa sample [1^1 , the BAO data from the SDSS Data Release 7 (DR7) galaxy sample [I3I , 
the CMB measurements given by the 7-year WMAP observations [5,], and the independent Hq measurements from 
HST [1^ (see Appendix A for details). With these observational data, the cosmological constraints on the ACDM 
model are obtained: n^ = 0.275l!]:J]^^ and h = 0.704t°:jj^^ at la confidence level (CL), while n^ = 0.275to:o38 and 
h = 0.70'ito°Qll at 2cr CL. 

Now we discuss the cosmic age problem associated with those 9 extremely old GCs in M31 galaxy. The results 
are given in Fig|T] Our method is as follows. Firstly, since in a fiat universe the age problem is mainly constrained 
by observational data of {ilm, h), by using the current cosmological observations (SNIa-|-BAO+CMB+i?o), we plot 
the la and the 2a confidence regions of the ACDM model in the ftm-h plane. Secondly, since a specific cosmic age 
corresponds to a specific diagonal line in the flm-h plane and a larger cosmic age corresponds to a lower diagonal line, 
by searching the diagonal line that tangents to the 2a confidence region of the ACDM model, we obtain the theoretical 
upper limit of present cosmic age (at 2a CL) predicted by the ACDM model. As seen in the figure, this theoretical 
upper limit is to < 13.88 Gyr. Notice that this upper limit is only given by the combined SNIa-|-BAO-|-CMB-|-7fo 
data, and has nothing to do with the age data of those GCs. Thirdly, we compare this theoretical upper limit of 
present cosmic age with the observational lower limits of those old GCs' ages listed in Table [TJ and find that the la 
lower limits of 5 GCs' age (including B024 with a lower limit 14.50 Gyr, B050 with a lower limit 15.70 Gyr, B129 
with a lower limit 14.40 Gyr, B297D with a lower limit 14.33 Gyr, and B495 with a lower limit 13.99 Gyr) are even 
larger than 13.88 Gyr. So the existence of these 5 GCs are in tension (over 2 a CL) with the current cosmological 
observations. This means that if the age estimates of these 5 GCs [3,[3] are correct, the cosmic age puzzle still remains 
in the standard cosmology. It should be mentioned that during this analysis, only 9 most discrepant GCs are used to 
quote statistical discrepancy. One may wonder what the discrepancy is when all objects are taken into account. By 
treating all those 139 GCs in M31 galaxy as data points and utilizing the x^ statistics, this issue is also checked. We 
find that the ACDM model can accommodate 127 GCs, and the statistical discrepancy disappears when all objects 
are included. 

Next we turn to the high-z cosmic age problem associated with quasar APM 08279 + 5255 at z = 3.91. The age 
of this quasar had been estimated by studying its chemical evolution. Based on the evolution of Fe/0 ratio from 
the X-ray observation, Ref. 9] gave a rough estimate igso = (2.0 ~ 3.0) Gyr at z = 3.91. Soon after, by using a 



ACDM Model 




0.24 0.26 0.28 0.30 0.32 0.34 



Q 



FIG. 2: High-z cosmic age problem at z = 3.91 in the ACDM model. The dotted elliptical line, the dashed elliptical line, and 
the central star symbol denote the la confidence region, the 2a confidence region, and the best-fit point of the ACDM model 
in the ilm-h plane, respectively. The solid diagonal line represents the lower limit of the quasar's age (i.e. 1.8 Gyr) at z = 3.91. 
The area below this diagonal line corresponds to a cosmic age larger than 1.8 Gyr. The ACDM model predicts a lower cosmic 
age than the la lower limit of the old quasar's age. 

detailed chemodynamical model for the evolution of spheroid, Friaca et al. [lOl obtained tqso = 2.1 ± 0.3 Gyr at 
the same redshift. The high- 2; cosmic age problem at z = 3.91 in the ACDM model is demonstrated in FigUl The 
method of analysis is same as that used in FigH] Notice that the solid diagonal line represents the la lower limit 
of the quasar's age (i.e. 1.8 Gyr) a,t z — 3.91. Our calculations show that the existence of quasar 08279 + 5255 is 
evidently inconsistent with the current cosmological observations, and the discrepancy is over 7a CL. 



III. COSMIC AGE PROBLEM IN THE INTERACTING DE MODELS 



As seen in the previous section, the ACDM model has difficulty to pass the cosmic age test. One may ask how other 
DE models fare. In a previous paper 1J|, Wang and Zhang demonstrated that by simply introducing the DE alone 
we cannot remove the high-^; age problem, and suggested that the introduction of interaction between dark sectors 
may be helpful to alleviate the cosmic age problem. So in this work, we also investigate the cosmic age problem 
in the interacting DE models. As examples, three kinds of interacting DE models are considered here. In the first 
subsection, we introduce these 3 interacting DE models, and determine their model parameters by using cosmological 
observations. Then in the second subsection, we evaluate the age of the universe in these 3 models, and use the 
obtained results to analyze the consistency with the age data listed above. 



A. Cosmological Constraints on Three Interacting DE Models 

The universe is filled with three classes of major energy components: matter (including baryons and dark matter), 
dark energy, and radiation. After introducing interaction between dark sectors, the dynamical evolutions of these 3 
components satisfy 



Pn 


+ iHpra 


= Q 


Pd 


, -f 3i/(pd, 


'.+Pde) 


Pr 


+ 'IHpr = 






-Q 



(7) 
(8) 
(9) 



where the over dot denotes the derivative with respect to the cosmic time t, Q denotes the phenomenological interaction 
term, and pde is the pressure of DE given by the specific model. 
We consider a simple scalar DE model with Lagrangian [19| 



1 



a'^</.a^0-y(0). 



(10) 



Assuming the potential energy V{4>) be dominant, the energy density and the pressure of DE are 

Pde = -Pde — V. (11) 

This means that the equation of state of DE w = Pde/pde = — 1- Noting that a = jr^, H = ^, we have 



d dt da 



dz dt da dz H{1 + z) dt 



(12) 



Then Eqs.(l7]l9|) can be rewritten as 



{l + z)^-3p,n^-Q/H (13) 

il + z)^ = Q/H (14) 

{l + z)^-Ap,=0 (15) 

dz 

Ea.([T5|) has a general solution pr = Proi^ + z)^, while the solutions of Eqs. (|13|) and ((T4l) depend on the specific form 
ofQ. 

Owing to lack of the knowledge of microscopic origin of the interaction, here we just follow ,2ft] and consider three 
classes of interaction terms 

Qi = "iaHpjn (16) 

Qs = iaHpde (17) 

Qs = iaH{p„, + Pde) (18) 

where a is a model parameter denoting the strength of interaction. For simplicity, we will denote these 3 models by 
lACDMl, IACDM2, and IACDM3, respectively. Notice that all these interacting DE models have 3 model parameters, 
Q.m, a, and h. 

First we discuss the lACDMl model. Substituting Qi — SaHpm into Eq. ([T^ . we have 

Pm=P™o(l+^)'^'""^. (19) 

Substituting Eg. ([TO)) into Eg. p^ and using the initial condition pde{z = 0) = pdeo, we get 



Pde 



"^"°(i + .)3a-)-^ + p,,o. (20) 



1 — a 1 — a 

Combining Egs. p^ . (1^0]) . and ([5]), we obtain the reduced Hubble parameter of the lACDMl model 



,) = Jn,{l + z)4 + i^(l + ;,)3(l-a) + h-nr- -^) 

V 1 — a 1 — a 



E{z) = A/r!,(l + Z)4 + -^^(1 + z)3(l-") + (l-nr- -^^ . (21) 

V 1 — a 1 — a 

Based on the x^ statistics method and observational data, the model parameters of the lACDMl model are determined 
as n„, = 0.279to;°27' " = (-1.20 x W^^)tlil'^ll~-l, and h = O.TOStao^g at la CL, while rj,„ = 0.279J:o:°g, 

a = (-1.20 X 10-^);5;oi^io-3, and /i = O.TObtofsl at 2ct CL. 

We now turn to the IACDM2 model. Substituting Q2 — 3aHpde into Eg. dTi)) . we have 

Pde = Pdeoil + zf"- (22) 

Substituting Eg. p2)) into Eg. p^ and using the initial condition p„i{z = 0) = p„M, we get 

Pm = /Omo(l + 2) - -; (1 + ^) +-; (1 + ^) ■ (23) 

1 — a 1 — a 

Combining Egs. ip^ . (I^5|) . and ([2]), we obtain the reduced Hubble parameter of IACDM2 model 



E{z) = \\^ri\ + ^)4 + ""' ^ + '^"'■ (1 + ^)3 + ^ "'" "- (1 + zf-. (24) 

" 1 — a 1 — a 



The model parameters of the IACDM2 model are determined as f],„ = 0.276lo;o2g, a = (-1.43 x 10"^)+i[-5g^io-2, 

and h = 0.702to:o24 at Icr CL, while Vi^ = 0.276to:o42> " = (-1.43 x lQ^^)tl\Hw-^, ^^^ h = 0.702^0.038 at 2cr CL. 
Next we consider the IACDM3 model. Substituting Qs = 3aH{pm + Pde) into Eqs. (IT5)) and (ITil) . we have 

(1 + z)--p- = 3(1 -a)p,„ - 3a/9de (25) 

{I + z)^ ^ 3a{p,-n + Pde) (26) 

Utilizing the initial conditions pmiz = 0) = pmo and pde{z = 0) = pdeO^ the differential equations (j25p and (j26p are 
solved as 

p„, = Ci(l + z)'=^+C2(l + z)'=% p,, = Ci ^~^"~^^ (l + z)^^+g /~^"~^^ (l + z)^-, (27) 

6a 6a 



where 



and 



3 



fci = -(l + VT^4^), fc2= 2(1-Vl^4^)- (28) 



^ _ (3 - 3q; - k2)pm0 ~ 3apdeO ^ _ (3 - 3q; - fci)jOmO - 'iapdeo .2q\ 

fcl - fc2 ' fc2 - fci 

From the Friedman equation, we obtain the reduced Hubble parameter of the IACDM3 model 



E{z) = ^/a.(i + zy + Ci^^ii + z)'^^ + C2^^(i + zY- (30) 



where 



- (3 - 3a - fc2)f^m - 3a(l - ri„j - f7r) ^ [i - ia - ki)Q.jn - ia[l - ^rn - ^r) /oiN 

L-i ; ; ) '-'2 ; ; . {61} 

Parameters of the IACDM3 model are determined as f7„ = 0.279to:o27> " = (-8.98 x 10-'')+2:?2xio-'' and h = 
0.704t°:g}^ at Ict CL, while Vt^ = 0.279tJ]°^^, a = (-8.98 x I0''^)ttl2xw-'^ ^ and /i = 0.704tS]:[!io at 2cr CL. 

B. Age Problem of Three Interacting DE Models 

In this subsection, we evaluate the age of the universe in these 3 interacting DE models, and analyze the consistency 
with the age data listed above. Since in a flat universe the age problem is mainly constrained by the observational 
data of (rim, h), our results are given in the flm-h plane. The method of analysis is same as that used in section IL 

We plot Figl3] to demonstrate the cosmic age problem in the lACDMl model. The left panel of Figl3] represents 
the cosmic age problem at z = 0. From this panel, it is found that in the lACDMl model, the present cosmic age 
^0 < 14.02 Gyr at 2a CL. Although the lACDMl model can give an upper limit of the present cosmic age 0.14 Gyr 
larger than that given by the ACDM model, this cosmic age upper limit is still smaller than the Icr lower limits of 4 
GCs' age (including B024 with a lower limit 14.50 Gyr, B050 with a lower limit 15.70 Gyr, B129 with a lower limit 
14.40 Gyr, and B297D with a lower limit 14.33 Gyr). So the lACDMl model still cannot pass the cosmic age test of 
these 4 GCs. We also check what the discrepancy is when all 139 GCs are taken into account. It is found that the 
lACDMl model can accommodate 127 GCs, and the statistical discrepancy disappears when all objects arc included. 
Besides, the right panel of Fig [3] represents the high-z cosmic age problem at z — 3.91. Our calculations show that the 
existence of high-z quasar APM 08279 + 5255 cannot be accommodated in the lACDMl model, and the discrepancy 
is at 5a CL. 

The situation of the IACDM2 model is demonstrated in the Fig|4l The left panel of Fig|4] represents the cosmic 
age problem at z = 0. From this panel, it is found that in the IACDM2 model, the present cosmic age to < 14.03 
Gyr at 2a CL. Although the IACDM2 model can give an upper limit of the present cosmic age 0.15 Gyr larger than 
that given by the ACDM model, this cosmic age upper limit is still smaller than the la lower limits of 4 GCs' age 
(including B024, B050, B129, and B297D). We also check what the discrepancy is when all 139 GCs are taken into 
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FIG. 3: Cosmic age problem in lACDMl model. The dotted elliptical line, the dashed elliptical line, and the central star 
symbol denote the la confidence region, the 2a confidence region, and the best-fit point of the lACDMl model in the Qm-h 
plane, respectively. In the left panel, three solid diagonal lines, from top to bottom, represent the constraints of 13.88 Gyr, 
14.02 Gyr, and 14.2 Gyr at z = 0, respectively. When plotting these diagonal lines denoting cosmic age constarints, the best-fit 
value of interaction strength of lACDMl model a — —1.20 x 10"'' is adopted. Notice that the top diagonal line tangents to 
the Ict confidence region of the lACDMl model, while the middle diagonal line tangents to the 2a confidence region of the 
lACDMl model. In the right panel, the solid diagonal line represents the lower limit of the quasar's age at z = 3.91. The area 
below this diagonal line corresponds to a cosmic age larger than 1.8 Gyr. The lACDMl model predicts a lower cosmic age than 
the la lower limits of 4 GCs' age, and also predicts a lower cosmic age than the la lower limit of the old quasar's age. 
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FIG. 4: Cosmic age problem in the IACDM2 model. The dotted elliptical line, the dashed elliptical line, and the central star 
symbol denote the la confidence region, the 2cr confidence region, and the best-fit point of the IACDM2 model in the Q,m-h 
plane, respectively. In the left panel, three solid diagonal lines, from top to bottom, represent the constraints of 13.91 Gyr, 
14.03 Gyr, and 14.2 Gyr at z = 0, respectively. When plotting these diagonal lines denoting cosmic age constraints, the best-fit 
value of interaction strength of the IACDM2 model a = —1.43 x 10~^ is adopted. Notice that the top diagonal line tangents 
to the la confidence region of the IACDM2 model, while the middle diagonal line tangents to the 2a confidence region of the 
IACDM2 model. In the right panel, the solid diagonal line represents the lower limit of the quasar's age at z = 3.91. The area 
below this diagonal line corresponds to a cosmic age larger than 1.8 Gyr. The IACDM2 model predicts a lower cosmic age than 
the la lower limits of 4 GCs' age, and also predicts a lower cosmic age than the la lower limit of the old quasar's age. 



account. It is found that the IACDM2 model can accommodate 127 GCs, and the statistical discrepancy disappears 
when all objects are included. Besides, the right panel of Fig|3] represents the high-z cosmic age problem at z = 3.91. 
Our calculations show that the IACDM2 model cannot pass the cosmic age test of the quasar APM 08279 -I- 5255, 
and the discrepancy is at 6cr CL. 

The situation of the IACDM3 model is shown in the FiglS] The left panel of Fig[5] represents the cosmic age 
problem at z = 0. From this panel, it is found that in the IACDM3 model, the present cosmic age tg < 13.94 Gyr at 
2(7 CL. Although the IACDM3 model can give an upper limit of the present cosmic age larger than that given by the 
ACDM model, this cosmic age upper limit is still smaller than the la lower limits of 5 GCs' age. We also check what 
the discrepancy is when all 139 GCs are taken into account. It is found that the IACDM3 model can accommodate 
127 GCs, and the statistical discrepancy disappears when all objects arc included. Besides, the right panel of FigjS] 
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FIG. 5: Cosmic age problem in the IACDM3 model. The dotted elliptical line, the dashed elliptical line, and the central star 
symbol denote the la confidence region, the 2cr confidence region, and the best-fit point of the IACDM3 model in the Qm-h 
plane, respectively. In the left panel, three solid diagonal lines, from top to bottom, represent the constraints of 13.84 Gyr, 
13.94 Gyr, and 14.2 Gyr at 2; = 0, respectively. When plotting those diagonal lines denoting cosmic age constraints, the best-fit 
value of interaction strength of the IACDM3 model a = —8.98 x 10~^ is adopted. Notice that the top diagonal line tangents 
to the la confidence region of the IACDM3 model, while the middle diagonal line tangents to the 2a confidence region of the 
IACDM3 model. In the right panel, the solid diagonal line represents the lower limit of the quasar's age at z = 3.91. The area 
below this diagonal line corresponds to a cosmic age larger than 1.8 Gyr. The IACDM3 model predicts a lower cosmic age than 
the la lower limits of 5 GCs' age, and also predicts a lower cosmic age than the la lower limit of the old quasar's age. 

represents the high-2 cosmic age problem at 2 = 3.91. Our calculations show that the existence of high-z quasar APM 
08279 -I- 5255 cannot be accommodated in the IACDM3 model, and the discrepancy is at 6a CL. 

Therefore, although the introduction of the interaction between dark sectors can give a larger cosmic age, the 
interacting dark energy models still have difficulty to pass the cosmic age test. This conclusion is different from that 

of ill. 



IV. SUMMARY 

The cosmic age problem is a longstanding issue, and provides an important tool for constraining the expanding 
history of the universe [21| . In this work, we investigate the cosmic age problem associated with 9 extremely old 
globular clusters in M31 galaxy and 1 very old high-2 quasar APM 08279 -t- 5255 at 2 = 3.91. It should be pointed 
out that these 9 GCs have not been used to study the cosmic age problem in the previous literature. By evaluating 
the age of the universe in the ACDM model with the constraints from the SNIa, the BAO, the CMB, and the 
independent Hq measurements, we find that the existence of 5 globular clusters and 1 high-2 quasar are in tension 
(over 2(7 confidence level) with the current cosmological observations. Therefore, if the age estimates of these objects 
are correct, the cosmic age puzzle still remains in the standard cosmology. Moreover, we extend our investigations to 
the cases of the interacting dark energy models. By studying the cosmic age problem in the three classes of interacting 
models, it is found that although the introduction of the interaction between dark sectors can give a larger cosmic 
age, the interacting DE models still have difficulty to pass the cosmic age test. This conclusion is different from that 
of [lj| . In a latest paper [23] , Cui and Zhang argue that the high-z cosmic age problem can be greatly alleviated when 
the interaction and spatial curvature are both introduced in the holographic dark energy model [2J, [2J|. Therefore, 
the cosmic age problem still needs to be further investigated in the future work. 
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Appendix A: Utilizing Various Cosmological Observations to constrain DE Models 



First we start with the SNIa observations. We use the Constitution sample including 397 data that are given in 
terms of the distance modulus fJ.obs{zi) [16|. The theoretical distance modulus is defined as 



^J'th{z^) = 5logio DL{Zi) + ^IQ, 



(32) 



where fiQ = 42.38 — Slogj^Q h with h the Hubble constant Hq in units of 100 km/s/Mpc, and in a flat universe the 
Hubble- free luminosity distance Dl = HqcIl {cIl denotes the physical luminosity distance) is 



Dl{z) = {1 + z) 



dz' 



Eiz'; 



where 6 denotes the model parameters. The x^ for the SNIa data is 



x|ArW=E 



ilJ'obsiZi) - ljLth{zi]0)f 



(33) 



(34) 



where ^obs{zi) and Ui are the observed value and the corresponding la error of distance modulus for each supernova, 
respectively. Following Ref. 25| , the minimization with respect to /iQ can be made trivially by expanding the x^ of 
Eq. (|34|) with respect to /.to as 



where 



[i^obs{zi) ~ iith{zi\y.o = 0,61)] 2 



A{e) = Y. 



cft 



(35) 



(36) 



r,,n\ _ V^ ^J'obs{zi) — flthiZi] fJQ — 0, 



Cf 



(37) 



^-Ei 



(38) 



Evidently, Ea.([M|) has a minimum for ^q — B/C at 



x%AO) = A{e)- 



B{B? 
C 



(39) 



Since 
Mo- 



XsN,t 



Xsiv mim instead minimizing X57V we minimize X57V which is independent of the nuisance parameter 



Next we turn to the BAO observations. The spherical average gives us the following effective distance measure [26 

.. li/3 



Dv{z) 



il + zYDiiz) 



where Da{z) is the proper angular diameter distance 



DAiz) 



H{z) 



dz' 



l + zj, E{z') 



(40) 



(41) 



The BAO data from the spectroscopic SDSS DR7 galaxy sample [Hl give Dv{z = 0.35)/Dv(z = 0.2) = 1.736±0.065. 
Thus, the x^ for the BAO data is. 



2 _ i'Dv(z = Q.?,b)/Dv{z = 0.2) - 1.736 
^^^° " I 0065 



(42) 
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Next consider the CMB observations. Here we employ the "WMAP distance priors" given by the 7-year WMAP 
observations (5|]- This includes the "acoustic scale" I a, the "shift parameter" R, and the redshift of the decoupling 
epoch of photons z*. The acoustic scale I a describes the distance ratio DA{z^)/rs{z^), defined as 

lA^{l + z.)^^^, (43) 

rs(z*) 

where a factor of (1 + z^) arises because Z?^(z*) is the proper angular diameter distance, whereas rs{z-^,) is the comoving 
sound horizon at z*. The fitting formula of rs{z) is given by 

rs(z) = —= / , (44) 

where 51;, and Vt-^ are the present-day baryon and photon density parameters, respectively. In this paper, we adopt 
the best-fit values, rJb = 0.022765/1-^ and n^ = 2.469 x lO'^/i-^ (for Tc„^b = 2.725 K), given by the 7-year WMAP 
observations b\ . The fitting function of z, is proposed by Hu and Sugiyama [23] : 

z, = lOAd,[l + {)mi24,{nbh'^)-°-'^^\l+gi{VLr„h'^Y% (45) 

where 

_ 0.0783(f}bfe^)-°-^38 ^ 0.560 

^^~ 1 + 39.5(f^6/l2)0.763 ' 92- ^^ 21.1(f^f,/l2)l-81 " ^ ^> 

The shift parameter R is responsible for the distance ratio _D^(z»)/iJ^^(z»), given by [28j 



R{z,) = ^n„,H§{l + z,)Da{z,). (47) 

Following Ref. 5], we use the prescription for using the WMAP distance priors. Thus, the x^ ^oi the CMB data is 

x'cMB = (^f - xf ^)(C-i),,(xf - xf ), (48) 

where Xi — {Ia,R,z^,) is a vector, and (C^^)y is the inverse covariance matrix. The 7-year WMAP observations [5| 
give the maximum likelihood values: Ia{z*) = 302.09, R{z^) = 1.725, and z, — 1091.3. The inverse covariance matrix 
is also given in Ref. Q 



(c-i) = 



^ 2.305 29.698 -1.333 

29.698 6825.27 -113.180 | . (49) 

»^ -1.333 -113.180 3.414 



At last, we also use the prior on the present-day Hubble constant Ha = 74.2 ± 3.6km/s/Mpc [l^. In Ref. |18| . the 
authors obtain this measured value of Hq from the magnitude-redshift relation of 240 low-z type la supernovae at 
z < 0.1. It is remarkable that this Gaussian prior on Hq has also been used in the analysis of the WMAP 7- year 
observational data @- The x^ function for the Hubble constant is 

2 /;i-0.742y , , 

^^ = hio3^ ■ (^«) 
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